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ABSTRACT 


The unsteady flow phenomena resulting from the interaction of wakes and 
vortices with the free surface are of particular importance in naval 
hydrodynamics. Ship and submarine wakes produce a three-dimensional 
complex signature, comprised of a narrow dark band bordered by two bright 
lines in synthetic-aperture-radar (SAR) images. The dark band signifies the 
suppression of waves at the Bragg frequency as a consequence of the interaction 
between the free surface and the imposed vorticity. In the present investigation, 
the vorticity field is provided by a single tip vortex generated by an airfoU. The 
results, obtained with an LDV, have shown that the free surface redistributes part 
or all of the normal turbulent kinetic energy into streamwise and spanwise 
components. The tiu*bulent kinetic energy first deaeases sharply with increasing 
vertical distance from the vortex and then remains nearly constant within a thin 
layer below the 'roughened' free surface. The results explain the longevity of the 
structures and lend further credence to the simulation of near-surface structures 
via vortex- or contour-dynamics. 
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I. INTRODUCTION 


The unsteady flow phenomena resulting from the interaction of wakes 
and vortices with the free surface are of particular importance in naval 
hydrodynamics. Ship wakes produce a three-dimensional complex signature, 
comprised of a ivarrow dark band bordered by two bright lines in synthetic- 
aperture-radar (SAR) images. The dark band is the most prominent of all the 
signatures and is seen many kilometers downstream at all angles to the SAR 
azimuth direction even under severe weather conditions. It signifies the 
suppression of waves at the Bragg frequency as a consequence of various 
short-wave-damping phenomena such as turbulence, surface-active 
materials, and the redistribution of surface impurities. The two bright lines, 
on the other hand, each resembling a moonglade, manifest themselves only 
in light winds and signify the occurrence of a range of waves which happen to 
be near the Bragg wavelength, possibly as a consequence of the interaction 
between quasi-two-dimensional turbulent motions near the free surface and 
the restructuring and modulation of this interaction by wind. 

The foregoing strongly suggests that the three-dimensional turbulent 
flow field beneath the flee surface, the extent of the two-dimensionalization 
of turbulence with depth, and the reverse energy cascade process near the free 
siirface must be understood in as much detail as possible in order to gain 
some iiisight into the occurrence of surface signatures. One of the simplest 
possible flows relevant to the dynamical processes in vorticity/free-surface 
interaction which can be carefully studied in isolation, without complications 
and competing influences that normally occur in a fully turbulent ship wake, 
is the interaction of a single turbulent vortex with the free surface of an 
otherwise smooth uniform flow. The modulations of the flow field and 
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tiirbulence near the free surface are not expected to be similar to that obser\'ed 
for a streamwise vortex in or near a rigid-wall boundary layer (Harvey and 
Perry, 1971; Shabaka, Me! ‘ a, and Bradshaw, 1985). 

Extensive reviews of the interaction of a pair of heterostrophic vortices 
with a free surface are given by Sarpkaya (1986, 1992a, 1992b), Sarpkaya, 
Elnits^', and Leeker (1988) and Sarpkaya and Suthon (1991), wherein it is 
suggested that the free surface signatures (scars) exhibit features that obey the 
decay laws of two-dimensional turbulence. The reverse energy cascade 
process dominates the turbulence behavior near the free surface through the 
interaction and merging of homostrophic vortices. The energy and enstrophy 
inertial ranges coexist, with an upscale energy transfer and downscale 
enstrophy transfer in the same wavenumber interval. The quasi-two- 
dimensionalization of turbulence in various types of flows (Jacquin, Leuchter, 
and Geffroy, 1989; Hunt and Graham, 1978), in general, and near a free surface 
(Loewen, Ahlborn, and Filuk, 1986; Brumley, 1984; Dickey, Hartman, 
Hammond, and Hurst, 1984; Komori, Ueda, Ogino, and Mizushina, 1982; 
Rashidi, and Banerjee, 1988; Lam, and Banerjee, 1992; Anthony, 1992; Komori, 
Nagaosa, Murakami, Chiba, Ishii, and Kuwahara, 1993), in particular, has 
emerged as a fundamental phenomenon during the past decade. 

In view of the foregoing, an experimental investigation on the flow 
structure resulting from the interaction of a single tip vortex with a 
deformable free surface was undertaken (Neubert, 1992; Sarpkaya and 
Neubert, 1993). The present study is a continuation of this investigation in 
much greater detail as far as the velocity and turbulence measurements are 
concerned. The results have shown that the free surface redistributes part or 
all of the normal turbulent kinetic energy into streamwise and spanwise 
components. The turbulent kinetic energy first decreases sharply with 
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increasing vertical distance from the vortex and then remains nearly constant 
within a thin layer below the 'roughened' free surface. The results lend 
further credence to the simulation of near-surface structures by a decaying 2- 
D turbulence through the use of vortex element methods. 
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II. EXPERIMENTAL EQUIPMENT 


The experiments were conducted in a low turbulence water turmel 
with an open test section of about 40 cm wide, 50 cm deep (maximum), and 
125 cm long. The turbulence management system is located upstream of the 
test section. It consists of a honeycomb and a fine-mesh screen. The tunnel is 
driven by a 50 Hp centrifugal pump. A second but smaller pump 
continuously circulated the tunnel water through a micro filtration system to 
remove rust and other suspended fine particles, down to about 1 pm, from 
the water (the filtration system was turned off during the experiments). A 
skimmer, placed upstream of the turbulence management system and 
coimected to another filter, was operated continuously to maintain the free 
surface clean. 

A 6% thick, rectangular planform, symmetric, Joukowski half-span foil 
with an effective aspect ratio of 10.4 was used to generate a 'single' trailing 
vortex. The tip of the foil was carefully rounded. The model foil was 
mounted in a rotatable cylindriced base, embedded into the bottom of the test- 
section floor. The angle of attack was varied from 6° to and the water 
level was held constant. The interior of the model was hollowed and 
connected to a dye reservoir to seed the vortex core with a fluorescent dye. 
The leading edge of the foil was 3.6 chord lengths downstream of the test 
section entrance. 

The mean velocities and turbulence intensities were measured with a 
Laser-Doppler Anemometer. Bragg-cell frequency shifting by 0.5 MHz was 
used in one channel to detect the flow reversals. The probe volume 
(approximately 200 pm in diameter) was positioned at the required location 
by use of a remotely driven x-y-z traversing unit. The scattering particles used 
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were titaiuum dioxide of rutile crystallme form and were approximately 12 
pm in size. All measurements were made using 1024 point e]:\sembles. The 
data are reported without any velocity-bias-correction. 

A single vortex position relative to the free surface and two 
measurement stations relative to the trailing edge of the foil were considered. 
The vortex axis was ho = 25 mm below the free surface. The velocity and 
hirbulence measurements were made along several y/ho = constant (vertical 
lines) and z/ho = constant (horizontal lines) at two stations along the vortex 
(5.0 and 6.5 chord lengths downstream from the trailing edge of the foil). The 
chord-based Reynolds number was UqC/v = 45,000, (F/v = 2,500), and the 
vortex circulation was about F/UqC = 0.50. Here only the most representative 
results will be discussed in terms of velocities and the RMS values of 
turbulence, denoted by u', v', and w* and normalized by Uq- 
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III. DISCUSSION OF RESULTS 

A. DISCUSSION OF VELOCITIES 

The coordinated axes and velocities are shown in Fig. 1. The axial 
velocity, normalized by the mean velocity Uq, is shown in Figs. 2 and 3 for 
z/ho = 0.58 and z/h© = 0.0, respectively, for x/c = 5.02. The axial velocity along 
representative vertical lines is shown in Figs. 4 and 5 as a function of z/ho for 
y/ho = -0.27 and y/ho = 0.0, for x/c = 5.02. Figures 2-5 show that the presence 
of the vortex leads to an axial velocity defect directly above the vortex axis. 
The axial velocity recovers quickly both along the horizontal and vertical 
lines as the values of either y/ho or z/ho increase. 

The vertical component of velocity, again normalized by Uq, is shown 
in Figs. 6 and 7 (for z/ho = 0.38 and z/ho = 0.0) and in Figs. 8 and 9 (for y/ho = 
-0.40 and y/ho = -0.27), all for x/c = 5.02. It is clear that the magnitude of the 
vertical component of velocity exhibits the usual vortex behavior near the 
vortex axis. This is particularly evident along the horizontal lines. Along the 
z axis, the magrutude of W/Uo decreases sharply as one approaches the free 
surface (Figs. 8 and 9). A representative plot of the lateral component of 
velocity V/Uo is shown in Fig. 10 for y/ho = 0.0. It exhibits the usual vortex 
velocity distribution along a vertical line. 

The axial velocity, normalized by the mean velocity Uq, is shown in 
Figs. 11 and 12 for z/ho = 0.85 and z/ho = 0.0, respectively, for x/c = 6.52. The 
axial velocity along representative vertical lines is shown in Figs. 13 through 
15 a function of z/ho for y/ho = -0.46, y/ho = 0.25, and y/ho = 0.65, for x/c = 6.52. 
Figures 11-15 show that the presence of the vortex leads to an axial velocity 
defect at or near the vortex axis. The axial velocity recovers quickly both 
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along the horizontal and vertical lines as the values of either y/h© or z/h© 
increase. 

The vertical component of velocity, again normalized by Uq, is shown 
in Figs. 16 through 18 (for z/h© = 0.85, 0.62 and z/ho = 0.0) and in Figs. 19 and 
20 (for y/ho = -0.45 and y/ho = -0.25), all for x/c = 6.52. Once again, it is clear 
that the magiutude of the vertical component of velocity exhibits the usual 
vortex behavior near the vortex axis. This is particularly evident along the 
horizontal lines. Along the z axis, the magnitude of W/Uq decreases sharply 
as one approaches the free surface. 

B. DISCUSSION OF TURBULENCE MEASUREMENTS 

The turbulence intensity u'/Uq = Vu^/Uq = Root-Mean-Square value 
of the u'-component of turbulence normalized by Uq is shown as a function 
of y/ho in Figs. 21 and 22 for z/ho = 0.60 and z/ho = 0.0, and in Figs. 24 and 25 
for y/ho = -0.27 and y/ho = 0.0, all for x/c = 5.02. It is dear from the last two 
figures that u'/Uo first decreases as z/hoincreases and then, very near the free 
surface, it increases sharply. Figure 25 shows v'/Uo as a function of z/ho ^or 
y/ho = 0.0 and x/c = 5.02. Once again, v’/Uq also decresaes with increasing 
2 /ho as one approaches the free surface. However, very near the free surface 
v'/Uo increases sharply, almost doubling its value. On the basis of these 
observations, an increase in w’/Uo will require an explanation for the 
increase in turbulent kinetic energy and a decrease in w'/Uo will signify the 
increase of anisotropy and the quasi-two-dimensionalization of an otherwise 
nearly isotropic turbulence. 

Figures 26 and 27 show representative horizontal profiles of w'/Uo at 
z/ho= 0.38 and at z/ho = 0.0 for x/c = 5.02. Figures 28 and 29 show 
representative profiles of w'/Uo at y/ho = -0.40 and at z/ho = -0.15 for x/c = 
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5.02. Clearly, w'/Uq decreases with increasing z/ho for a given y/ho- 
However, unlike u’/Uq and v‘/Uo, w’/Uq does not exhibit an increase near 
the free surface. This signals an interesting metamorphosis as far as the effect 
of the free surface on turbulence is concerned. The turbulence becomes much 
more anisotropic and the reduction in w'/Uq (the reflection of vertical 
momentum from the free surface) leads to an increase in u'/Uq and v'/Uq. It 
is, therefore, important to explore the variation of the turbulent kinetic 
energy, particularly in the region where the two-dimensionalization of 
tiurbulence occurs. 

Figures 30 and 31 show u’/Uq and w’/Uq versus z/ho for y/ho = -0.46 
and y/ho = -0.25, respectively, all for x/c = 6.52. Once again, it is clear that 
while u’/Uq begins to increase near the free surface, the w’/Uo component of 
turbulence continues to decrease. There is no doubt that v'/Uq exhibits a 
similar behavior at x/c = 6.52 as it did at x/c = 5.02. 

The foregoing shows that the vertical component of turbulence 
decreases in planes close to the free surface. Normally, one would assume 
that this should be so and that the vertical velocity and turbulence 
fluctuations must go to zero as a consequence of the free surface boundary. 
However, the interaction of the vortex with the free surface gives rise to 
considerable surface signatures (scars and whirls with vertical axes), and the 
vertical fluci nations need not go to zero as the surface is approached. 
Nevertheless, one would expect diminishing turbulence intensities, with 
finite terminal values, as the free surface is approached. 

Figure 32 shows the turbulent kinetic energy E = (u'2 + v'2 + w’2)/2, 
normalized by the minimum kinetic energy in each investigator's data, as a 
function of z/ho. As one might expect, E/Emin is near maximum on the 
surface-normal passing through the vortex center. As the free surface is 
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approached, E/Emm decreases sharply at first to about 0.018 (near z/ho = 0.6), 
then remains nearly coi\stant up to about z/ho = 0.9, and then shows a slight 
increase as z/ho-> 1 (the free surface). In other words, the two- 
dimensionalization of turbulence very near the free surface is not at the 
expense of total turbulent kinetic energy and that the increase of both u’/Uo 
and v’/Uo is indeed due to the reflection of vertical momentum from the free 
surface. In this reflection, the free surface behaves as a scarred surface, 
exhibiting hydrodynamic roughness as it is seen from below. Had it behaved 
like a rigid, shear-free boundary, the reflection process would have been 
relatively simple to interpret. However, the formation of a scar, on the 
down-flow side of the vortex, the unsteady fluctuations associated with the 
formation, amalgamation and motion of the whirls (small vortices with axes 
in the surface-normal direction), and, possibly, the interaction of these scar 
structures with the fluid (air) directly above render the reflection process 
much more complex than that in a fully-developed turbulent channel flow 
(Komori et al., 1982; Rashidi and Banerjee, 1988; Lam and Banerjee, 1992) and 
require a closer examination of the quasi-coherent turbulent structures near 
the free surface. 

The foregoing raises a number of additional questions which must be 
addressed if one is to understand the two-dimensionalization process and its 
role on the formation of SAR imdr,as: does E/Emin indeed increase very 
near the free surface, are there other flow situations (not involving a trailing 
vortex) which exhibit similar characteristics near the free surface, and how do 
these structures decay with time or downstream distance. Undoubtedly, these 
questions are rather complicated and their numerical resolution will pose 
severe tests for all turbulence models. 
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The transition from fully three-dimensional towards two-dimensional 
tiirbulence, the slowing down of the decay of kinetic energy and the strong 
increase of the length scale have attracted attention in other tiirbulent flows 
for a number of years. Jacquin et al. (1989) and the numerous references cited 
therein dealt with the effect of rotation on the two-dimensionalization of 
turbulence and the strong departure of the length scales from those for an 
isotropic state (Brumley, 1984). Hunt and Graham (1978) investigated 
theoretically the case of homogeneous free stream turbulence impinging 
upon a rigid wall and found a growing viscous boundary layer adjacent to the 
free surface and a larger inviscid 'source layer.’ At the edge of the viscous 
surface layer the kinetic energy of turbulence is the same as in the bulk of the 
fluid, but the fluid velocity component normal to the surface vanish and the 
energy in this motion is partitioned (equally, for an isotropic flow) to the 
streamwise and lateral motioi\s. 

Loewen et al. (1986) investigated the statistics of free-surface flow 
structures generated by a vertical bar grid moving through water in a towing 
taitic. They have found a profusion of coherent surface structures that either 
rotate (surface eddies), translate (river flow), or are relatively stagnant and 
that the surface flow is predominantly two dimensional. Their results show 
that the eddy size distribution gradually shifts to larger size as turbulence 
decays and an increasing fraction of the fluctuating kinetic energy is 
transferred into 'rivers'. Brumley (1984) and Dickey (1984), made turbulence 
measurements near the free surface in grid-stirred tanks and demonstrated 
that the true dissipation rate is relatively uniform near the surface, rather 
than decaying as the fourth p>ower of the distance from the horizontal grid, as 
it does in the case of deeply-submerged grids. Sarpkaya and Suthon (1991) 
suggested, on the basis of their observations, that the free surface turbulence 
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exhibits features that obey the decay laws of two-dimensional turbulence. 
Subsequently, Sarpkaya (1992a, 1992b) has simulated the evolution of near¬ 
surface structiires by a decaying 2-D turbulence through vortex dynamics and 
has shown that the energy and enstrophy inertial ranges coexist, with an 
upscale energy transfer (reverse energy cascade) and downscale enstrophy 
transfer in the same wavenumber interval while conserving energy and 
enstrophy. 

Komori et al. (1982) made measurements of temperature, velocity, and 
turbulence fields very close to the free surface in an open channel flow and 
found that the vertical motions are damped while the streamwise and 
spanwise motions are promoted. Subsequently, Rashidi and Banerjee (1988) 
and Lam and Banerjee (1992) found, through flow visualizations of bubble 
streaks near the free surface of an open-channel flow and also through 
numerical simulations, that the kinetic energy of turbulence is redistributed 
from the surface-normal to the parallel fluctuations. Their results are quite 
similar to those obtained by Komori et al. (1982) in the range where they may 
be compared. Lam and Banerjee (1992) noted that unlike the numerical 
prediction, the experimental values of w' did not vanish at the interface. 
They have suggested that "this was probably due to the slight waves present 
in the experiments." As noted earlier, these experiments are for relatively 
simple chaimel flows, containing no trailing vortices or jets. 

More recently, Anthony (1990) and Anthony and Willmarth (1992) 
made turbulence measurements in a round jet beneath a free surface and 
found that near the surface u'/Uq and v'/Uq increase while w'/Uq decreases. 
Their E/Emin data (Anthony and Willmarth, 1992), calculated from their Fig. 
4b, are shown in Fig. 32 together with those calculated from the RMS data of 
Lam and Banerjee (1992), Komori et al. (1982) and the present investigation. 
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In doing so, the vertical distance (from the channel bottom, the jet axis, or 
the vortex axis) was expressed in terms of the relative distance z/ho. The 
turbulent kinetic energy for each case was scaled by its corresponding 
minimum value. This normalization does not imply identical scaling laws 
but allows one to compare the decay of E/Emin for the four sets of data 
considered herein. 

The conclusion common to all of the foregoing studies is that the free 
surface redistributes part or all of the normal hubulent kinetic energy w'2 /2 
into the streamwise and spanwise components u’2/2 and v'2/2. The 
tiubulent kinetic energy increases very near the free surface in all of the four 
sets of data. The results for the channel and vortex flows [Lam and Benarjee 
(1992), Komori et al. (1982,1993), and the present data] show for z/ho > 0 that 
the values of E/Emin decrease by comparable rates and then acquire nearly 
constant values in a very thin layer near the free surface. This lends credence 
to Sarpkaya's (1992a, 1992b) and Sarpkaya and Suthon’s (1991) numerical 
simulation of the evolution of near-surface structures by a decaying 2-D 
turbulence via vortex dynamics (see, e.g., Sarpkaya, 1989). The case of the 
round jet parallel to the free surface is an exception where the turbulent 
kinetic energy remains nearly constant throughout the region above the jet 
axis. Apparently, the dynamics of turbulence very close to the free surface (as 
z/ho-»l) is very complex and sigruficantly different from that near a smooth 
or rough rigid surface. The role of the near-surface layer, the persistence of 
u'/Uo and v'/Uo as z/ho-^1, the interaction between quasi-two-dimerwional 
coherent structures near the free surface and the restructuring and 
modulation of this interaction by wind and contaminants need extensive 
measurements. 
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IV. CONCLUSIONS 


1. The interaction of a single turbulent trailing vortex with a deformable 
free surface has been undertaken for the purpose of exploring the velocity and 
turbulence field in the vicinity of a clean free surface. The evidence presented 
herein shows that there are two regions wherein one, near the vortex, the 
turbulent kinetic energy decays rapidly with increasing vertical distance, and 
in the other, very near the free siurface (in the order of millimeters), the said 
energy remains essentially constant. The energy of the vertical velocity 
fluctuations is redistributed among the lateral and streamwise motions. Even 
though all four flows discussed herein exhibit, to varying degrees of intensity, 
the behaviors noted above, trailing vortices and open channel flows delineate 
the two distinct regions of turbulent energy distribution more emphatically 
than submerged jets. Finally, it is noted that the results of the present 
investigation lend further credence to the numerical simulation of the free 
surface structures via two-dimensional vortex dyi\amics. 

2. The use of a turbulent vortex near the free surface of an otherwise 
smooth uniform flow proved to be a 'kernel' experiment towards the 
elucidation of the dynamical processes in vorticity/free-surface interaction. 
Detailed measurements of mean velocity and all components of the Reynolds 
stress tensor are needed towards the development of a predictive model of 
the vorticity/free-surface interaction. These measurements and the 
imderstanding of the character of the quasi-coherent structures near the free 
surface constitute the essence of the on-going investigation. 
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Figure 1. Definition Sketch 
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Figure 3. U/Uo versus Lateral Distance for Z/ho = 0.00 and x/c = 5.02 
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Figure 4. U/Uo versus Vertical 
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Figure 5. U/Uo versus Vertical Distance for Y/ho = 0.00 and x/c = 5.02 
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Figure 7. W/Uo versus Lateral Distance for Z/ho = 0.00 and x/c = 5.02 
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Figures. W/Uo versus Vertical 


























oil 



24 


x/c = 6.52 
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Figure 12, U/Uo versus Lateral Distance for Z/ho = 0.00 and x/c = 6.52 
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Figure 13. U/Uo versus Vertical 
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Figure 15. U/Uo versus Vertical 
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Figure 16 . W/Uo versus Lateral Distance for Z/ho = 0.85 and x/c = 6.52 
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Figure 17. W/Uo versus Lateral Distance for Z/ho = 0.62 and x/c = 6.52 












Figure 18. W/Uo versus Lateral Distance for Z/ho = 0.00 and x/c = 6.52 











Figure 19. W/Uo versus Vertical Distance for Y/ho = -0.44 and x/c = 6.52 












Figure 20. W/Uo versus Vertical Distance for Y/ho = -0.25 and x/c = 6.52 
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Figure 23. u'/Uo versus Vertical Distance for Y/ho = -0.27 and x/c = 5.02 
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Figure 24, u7Uo versus Vertical Dis 
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Figure 26. w'/Uo versus 
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Figure 27. w'/Uo versus 













Figure 28. w /Uo versus Vertical 
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Figure 29. w'/Uo 
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Figure 30. u7 Uo versus Vertical 
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Figure 31. w'/Uo 











Figure 32. Normalized Turbulent Kinetic Energy versus Vertical Distance 
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